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Abstract 

We calculate rates for hybrid meson production by electromagnetic and weak interactions 
in the flux-tube model. Applications include photo and electroproduction at Jefferson Lab- 
oratory and HERA, and the production of light strange and charmed hybrids in the weak 
decays of heavy flavours. Photoproduction of some light hybrids is predicted to be prominent 
in charge exchange reactions, 7p nTi. and accessible in yp — > pTl. Production of light or 
charmed hybrids in B and D decays may be feasible with high statistics. Photoproduction 
of the sxial hybrid meson is predicted to be large courtesy of vr exchange, and its strange 
counterpart is predicted in B ^ iPKh{1^) with b.r. ~ 10~*. Production rates for exotic 
hybrid candidates 1""*"; (0, 2)^~ are given special attention. Selection rules that can help to 
distinguish between hybrid and conventional states with the same J^'~" are noted. 
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Introduction 



An outstanding problem in the Standard Model is how the non-Abelian, gluon, degrees of freedom 
behave in the limit of strong QCD. Lattice QCD predicts a spectroscopy of glueballs ^ and hybrid 
mesons j2|, but there are no unambiguous signals against which these predictions can be tested. 

A major stumbling block in the case of hybrids is that while predictions for their masses |2IISI, 
hadronic widths |4I5| and decay channels |4l5l6j are rather well agreed upon, the literature contains 
no general discussion of their production rates in electromagnetic or weak interactions. Meanwhile 
a significant plank in the proposed upgrade of Jefferson Laboratory is its assumed ability to expose 
the predicted hybrid mesons in photo and electroproduction. 

Clearly a calculation of hybrid photoproduction in a model based on lattice QCD is urgently called 
for. Theory jJICl has provided compelling arguments from QCD that confinement occurs via the 
formation of a flux tube. In the simplest situation of a long tube with fixed Q, Q sources on its 
ends, a flux-tube has a simple vibrational spectrum corresponding to the excitation of transverse 
phonons in its string-like structure. There is a question whether the flux-tube is fully formed on 
the 0(1) fm scale typical of hadrons (ref. 8 suggests that at 0(1) fm the state of confined gluon 
fields is between that typical of bag models and those of a fully formed fiux tube; by contrast |5] 
finds that the flux-tube forms at distances below 1 fm) . The model of ref. 3 assumes that the fully 
formed tube drives the phenomenology and that the essential features of this gluonic spectrum 
are retained in the spectrum of real mesons with their flux-tube excited - the hybrid mesons. 

The flux tube model has not been derived from QCD, but it has many features that appear to 
be shared by lattice QCD. As such it is currently the best that we have that builds on features 
emerging from the strong coupling limit of QCD. Hopefully its successes or failings can enable 
deeper insights to emerge as to the dynamics of QCD in the strong interaction limit. It leads to 
the effective linear potential of the conventional meson spectroscopy and has an identical N=l 
spectrum of hybrid states as the lattice - a feature that is especially relevant to the present 
calculations. Its pedigree includes predictions that have inspired studies on the lattice and even 
anticipated their results (such as the J^'^ pattern of the lightest hybrids P], and that the decays 
of those hybrids to ground state mesons are suppressed relative to those to excited states |3] [5], 
which has recently been confirmed in the lattice framework ^Ol)- We would hope that insights 
from the present work in their turn might inspire future lattice study, which could thereby better 
establish their connection with reality. 

The model is widely used in the experimental community for whom the reasons underlying the lat- 
tice results are often obscure. As such it is playing a prominent role in the the Jefferson Laboratory 
upgrade proposal, where the flux-tube idea has been used to underpin much of their planning. To 
turn this into something of practical use will require predictions for the electromagnetic transition 
amplitudes to hybrid mesons. However, while implications for spectroscopy and hadronic decays 
in such a model have been extensively explored, previous estimates of electromagnetic couplings 
in this model (refs. ^2E1) are at best only upper limits, in that they were based upon vector 
meson dominance of hadronic decays of hybrids into modes including p and further assumed that 
the predicted suppression of decay into -np is suspended in tt exchange. We are unaware of any 
direct calculation of photo- or electro-production of hybrids in this model. This is the issue that 
we address in this paper. 

From among the results of our extensive survey we note the following. 

(1) The electric dipole transitions of the hybrid axial meson to 7r*7 and of the exotic (0, 2)^ to 
P7 give radiative widths that can exceed 1 MeV. This implies significant photoproduction rates 
in charge exchange reactions 7p — > H~^n. The exotic 2"' may also be produced diffractively in 
7P Hp. 

(2) The production of an axial strange hybrid in B ^ i^Kn is predicted to have b.r. ^ O(10~^) 
so long as its mass < 2.1GeV. There is a tantalising unexplained enhancement in the data that 
may be compatible with this (IS) and merits further investigation. 

(3) If there is a light exotic hybrid | 14II15[I16[[T71 with J^^ = 1"+, the 1=1 and 1=0 states could 
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be produced in Z? — )■ ttH and Ds — » vriJ with branching ratios ^10 '. 

(4) Measurement of the axial to vector amphtude ratio in B ^ D(^h)^ nia-Y enable the hybrid 
content of charmed mesons to be determined. 



1 The model 

The flux-tube is a relativistic object with an infinite number of degrees of freedom. A standard 
approximation p^Ell^lT^ has been to fix the longitudinal separation of the QQ = r and to solve 
the flux-tube dynamics in the limit of a thin relativistic string with purely transverse degrees of 
freedom. The resulting energies E{r) are then used as adiabatic effective potentials on which the 
meson spectroscopies are built. Ref. 19 studied the effect of relaxing these strict approximations 
and found that the spectrum of the conventional and lowest hybrids is robust. We shall assume 
the same is true in this first calculation of electromagnetic excitation of hybrid mesons. 

In refs. O^IE] the flux-tube was discretised into + 1 cells, (modern lattice computations 
typically have N ^ 10) and then N oo. Up to N modes may be excited. We shall focus on 
the first excited state, with excitation energy — where r is the length of the flux tube. If the 
length of a cell is I then r — {N + 1)1. 

The state of the fiux tube can be written in terms of a complete set of transverse eigenstates 

\y) = \yi .■■yn---yN) 

and the Fourier mode for the first excited state^ is 



- N 



sm (1) 



AT + 1 zl^ — iV + 1 

n=0 



or 



Vn — \ a sm (2 ) 

The oscillations are in the two dimensional space transverse to the nominal QQ axis. Thus there 
are two Fourier modes a = (ai, 02) where 1, 2 refer to the two (body-fixed) orthogonal coordinate 
directions, 61,62- 

In the small oscillation approximation the system becomes harmonic in y (a). Then if b is the 
string tension (6 1 GeV/fm), the eigenfunctions for the ground and first excited states (labeled 
0, 1 respectively) are in the Fourier-mode space 



1/4 

Xo(ai.2) = ( ] exp 



bn 2 



2{N + 1) 



(3) 



Xi(ai,2) = \J ^^^i °^i.2Xo(qi.2) (4) 

To reduce the number of indices write a = (01,02) in the body- fixed basis and understand any a 
to refer always to these components. To proceed to the continuum, write (see 3, 18 ) 

Pi 



N + 1 



The Gaussian wavefunctions eqns.(|31 Q for the flux-tube ground-state and first excited mode 
become 



^Higher modes (up to p = N) exist but we are not interested in them here as we wish to speciaUse to the hghtest 
hybrid mesons. Incorporating such modes is straightforward in the formahsm we describe. 
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Xi(a) = V2/3iaxo(a) (6) 

The wavef unctions for mesons must include the state of the flux-tube. So for conventional mesons, 
where the flux-tube is in its ground state, write 

C = V'ih„(^)Xo(ai)xo(a2), 

where n,l,m are the usual two-body quantum numbers and the subscript zero indicates the ground 
state. If either of the transverse modes is excited, one has a state that we refer to as a hybrid 
meson. The particular combinations '^[O'l i 202] = "Tf*^^ ^^^^ normalised circularly polarised 
phonon modes for the flux-tube, which have angular momentum ±1 about the longitudinal axis. 
The corresponding wavefunction for such a hybrid may be summarised by 

^ = V-itrnW-^ [Xi(ai)xo(a2) ± «Xo(ai)xi(a2)] 



or 



^ = Pi (ai ± ja2)xo(ai)xo(a2)- 



The challenge is : how can electromagnetic or weak currents, which couple to quarks, break the 
orthogonality of (a) needed to give a transition between conventional and (first) excited flux-tube 
(hybrid) state? The answer is implicit in the observation of Isgur that the flux tube is a 
dynamic entity, with a zero-point motion which can affect observables that at first sight are driven 
by the Q or Q. 

The physical picture becomes transparent if one simulates the tube as a series of beads with mass 
m on a massless string, and simplifies to the cases of iV = 1,2 [50]. We will be interested here in 
excitation of the first hybrid mode, whose essential spatial structure can be simulated by a single 
bead of mass m = br. 

In this first excited mode the centre of mass of the QQ-bead system is displaced from the interquark 
axis by a transverse distance that scales as m/niQ. If the transverse displacement of the bead 
is y, and the Q and Q have masses ttiq, then relative to the centre of mass, the position vector of 
the quark (antiquark) has components in the longitudinal r and transverse y directions 



'QiQ) 



f br 



"2 ' V 2to 



(7) 



The dependence of ?'q(q) on y enables a quark-current interaction to excite transitions in the y 
oscillator, leading to excitation of the fiux-tube. The presence of ± in the r coordinate, but only 
+ in the case of y(a), is a feature of the first excited mode. This is illustrated for the full flux-tube 
in flg.l^. The longitudinal axis passes through the cm. of the system. If the tube's effective cm. 
(the "bead") is displaced transverse to this in one direction, then the Q and Q respond collectively 
to the displacement of the flux tube, and are both on the opposite side of the longitudinal axis. 
Hence the same sign appears in the y{a) coordinate, but opposite signs in the longitudinal ±r. 
This sign will have signiflcant consequences when we discuss El transitions. 

This is the essential physics behind the excitation of hybrid modes by current interactions with 
the quark or antiquark. 

Extending to N beads [T^ leads to more mathematical detail and enables excitations of up to N 
modes, but the underlying principles are the same. The position vector becomes 



1^ ^ br 
I 



rQ^^^^R± r + d—-^l^^ (8) 
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Figure 1: p = 1 mode hybrid structure. rriT 



rrid 



br. 



or in the continuum hmit 



^Q(Q) = ^ ± 




(9) 



The essential physics is aheady contained in the above examples. For modes with p = even, the 
moment of the tube deformation relative to the interquark axis tends to cancel, leading to a null 
displacement of the quarks. For p — odd, the tube has a net transverse moment, leading to a 
compensating transverse displacement of the quark(s). This is encoded in the factor (— 1)^ for the 
quark displacement in the p-th mode in Isgur's formulation (eqn.(16) of _18 ). 

Isgur |18| noted that the dynamical degree of freedom implicit in the y (or a) gives a hitherto over- 
looked ( but welcome) contribution to the charge radius in elastic form factors. This comes about 
because the Coulomb interaction with the quark Fci{q) — {g.s.\e'^^'^'i\g.s.) receives contributions 
from both the r and the a degrees of freedom at 0{q^r^) and Oiq^a^). 

In the contimmm limit, N oo this becomes (compare eqns.(29) et seq in Isgur) 



k1 



24 



(1 



1 



(10) 



where the sum is over the p = 1 • • • oo "phonon" modes (elsewhere in most of this paper we 
consider only the p — I mode). It was shown JHI that these "transverse excursions" give large 
~ 51% corrections in light quark systems where mg = rud, and ~ 13% corrections in heavy-light 
Qq systems. Furthermore the J^Ti^/P^) ^ ^^"^ saturated by its p = 1 term. Together, these 



suggested that the transition amplitudes to the lowest hybrids (p 
substantial. 



1 phonon modes) could be 



This is our point of departure, 
momentum transfer, 



Expanding the incoming plane wave to leading order in the 



exp 



-iq ■ 





iq - a Pi 




(11) 



the linear terms in f and a break the orthogonality of initial and final wavefunctions and cause 
transitions: El excitation to L = 1 conventional states in the case of f and excitation of the first 
flux-tube mode (hybrid) in the case of a. By combining the above with the tensor decomposition of 
the current-quark interaction, we may calculate electromagnetic and weak excitation amplitudes 
to hybrids, and compare with those for conventional mesons in various multipoles. 

First we take the continuum limit, define wavefunctions for both hybrid and conventional states 
including the flux-tube, illustrate the familiar El electromagnetic transition and then calculate its 
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analogue for hybrid excitation with specific reference to the angular integrations. The electroweak 
transitions of heavy flavours will then be described, this requires knowledge of the decomposition 
of Vfj. and for heavy light systems. 

When niQ ^ mQ Isgur's decomposition of the position vectors was 



R±r- 



fl2) 



where ^ = mQirid/ [mq + to^). 

For nomenclature we shall adopt the PDG notation, where the subscript J denotes the total 
angular momentum of the meson, and we append the subscript H to denote hybrid (thus for light 
flavours we have the notation in Table. J^l, with obvious generalisation for flavoured states). This 
is done for two reasons 

(i) to enable the trivial distinction between hybrid and conventional states to be immediately 
apparent and reduce confusion in the text; 

(ii) as a reminder that for a conventional and hybrid meson with the same overall J^'^ ^ their 
internal qq spin states are inverted. 



IPC 


Sqq 


/ 1 


/ = 


1++ 





aiH 


flH 







PH 


Ll>h 


(0,1,2)+- 


1 


bjH 


hjH 


(0,1,2)-+ 


1 


T^JH 


VJH 



Table 1: Naming convention for light quark hybrid mesons 



This spin inversion is also illustrated in Table. (|T]l and has potentially important implications in 
helping to isolate hybrid contributions to the wavefunction of conventional states. For example, 
whereas a meson with J^'^ — 1 '";(0,2)"' has exotic correlations of J^'~' inaccessible to conven- 
tional mesons, the other states could a priori be conventional or hybrid. Note however that the 
1^^ have Sqq = whereas their conventional counterparts 1++ have Sqq = 1. Conversely, the 
(0, 2)]^+ and iff all have Sqq = 1 whereas their conventional counterparts all have Sqq = 0. Hence 
there is a complete spin inversion between conventional states and their hybrid counterparts. 

This spin inversion enables a dynamical distinction to ensue between these two types of state, which 
is manifested in certain selection rules. This has already been noted in hadronic decays and will 
have consequences in current transitions too. Furthermore, we will find that the hierarchy of spin 
operators leading to certain heavy flavour decays is inverted between conventional and hybrids, 
such that observable consequence can ensue in principle. 



2 El Transitions 

2.1 Conventional QQ Transitions 

To establish notation and make subsequent analysis of hybrid excitation more transparent, we 
first illustrate conventional El transition in the non-relativistic limit. Consider the El transition 
operator 

q=Q.d 

whose essential structure when acting on meson Qd is 
Oei ^-i\q] 
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U l,m denote the orbital angular momentum of the Qd system, and its ^;-projection (on fixed 
space axes) respectively, then the structure of the matrix element A4 becomes 

M ^ {l',m'\OEi\l,m) = -z\q] Srj d^a (r)x5(ai)xS(«2) e+-f (f)xo(ai)xo(a2) 

The integration over d?a ^ \ and the standard integral over (fir gives a transition from Z = to 

I' = \ caused by the presence of r. 

Separating into radial and angular parts, tpnimi'f^ = Rn{,r)Y"'- (Q) and noting that 



M. becomes 



f ^ - ^) fir). [ dn Yp'*Y+'Y{^ 

\mQ mdj V 3 J 



where f{r)i = J r^ dr R'^, (r) r Rn{r) . 

This general formula becomes more transparent when applied to the case Z = 0, T = 1 for which 

M = -i\q\^i - ^) f{r)~Sm',+i 

\mQ mdJ V3 

We are interested in the specific case: 

^bi) = i b{r),,\qb\^Sm',+i- (13) 

In general we can write the radiative width as 

r{A - Bj) = 4^|g1— i— ^ \M{mj,mS = mj + l)f (14) 
ruA 2,Ja + 1 

where the sum is over all possible helicities of the initial meson, and the matrix element is under- 
stood to be for a positive helicity photon. 

This corresponds to the familiar El transition formalism of atomic and nuclear physics as tradi- 
tionally applied to QQ systems. 



Notice that ( — — ) ensures charge-conjugation conservation; for charge-neutral systems the 

QQ charges cancel but they are vectorially on opposite sides of the cm. ("longitudinal" electric 
dipole moment). Hence a non- vanishing El amplitude occurs between neutral systems (e.g. x ~^ 
7V')- 

2.2 Transitions to Hybrids 

Transitions to hybrids in the first excited state of the flux-tube arise from the a component of the 
Q, Q position operators. 

We consider the general matrix element for transitions between a conventional meson and a first 
excited mode hybrid with tube oscillation polarisation ±1; 

M = (hyb; ±, m'|0|conv; l,m) = J (f f j d^a H*OC 

where (i) O - the essential spatial structure of the transition operator, which in this example is 



Oei = -i|<7l [^ + —]\^0ire+-a (15) 
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(ii) C - the conventional meson wavcfunction, 



<^ = ^itojW Xo(ai)xo(a2) 
(iii) TL* - the hybrid wavefunction (complex conjugate), 

^* = V'itl*(^) /^i («i T ia2) xS(«i)xS(a2) 

(where the flux tube is excited into a state with polarisation ±1 along its axis). 

To be specific, we consider the transition between the unexcited tube and a tube that has polari- 
sation ± along its axis (the "body axis"), that axis in turn being oriented at some angles 0,0 in 
the laboratory (the "fixed axes"). 

Later we will consider both vector and axial currents in various multipoles. In general, to have a 
transition between a normal meson and a hybrid, one power of a will be needed in the transition 
operator. The factors multiplying a will depend on the tensor structure of the current (e.g. vector 
or axial, transverse or longitudinal), explicit forms being given in Appendix^ Having identified 
the presence of a we need to be able to compute its expectation value. This is outlined in Appendix 
121 for a "reference" operator Oj-of = a ■ Xi. The El case then follows immediately. 

The integral over flux-tube variables gives 

(Xi,±|a.f±|xo) = ±j^(5+V^i': - (16) 
(Xi,±|a-z|xo) = -^(S+V'il* ^S-V'-^I*) (17) 



and so, for transition from a conventional ^5*0 state, for example, we have 



(hyb; ±, m'| exp 



-iq ■ r 



a ■ a;_ |conv; I — 0) 



1 



-^conv (^) 



(18) 



First expand the exponential in terms of partial wave angular states, contract together the three 
T) functions and integrate J c?f2, which gives for the matrix element (see Appendix Ib)| 



where we now only need to calculate the radial expectation values of the spherical Bessel functions. 

In the previous equations the factors 6^ refer to the flux tube polarisation transverse to the body 
vector r, while the 5m' ,±i refers to the meson's total angular momentum projection in the fixed axes 
(x, y, z) The parity eigenstates in the flux tube are given in ref. P). They are linear superpositions 
of states where the flux-tube has polarisation ±1. Following that reference we denote the number 
of positive or negative helicity phonon modes by which for our present purposes will be 

{1, 0} or {0, 1}. Parity eigenstates are then the linear superpositions 

|7'-±)^^(|{1,0})T|{0,1})) 

The effect is that when we take expectation values for parity eigenstates, the terms proportional 
to {5~^ ± 5" ) will be destroyed for the "wrong" parity, and amplified by \/2 for the "correct" parity. 
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This is the source of the extra overaU factor of \/2 in the foUowing expressions for transitions to 
specific parity eigenstates. With this preamble we now proceed to complete the expression for 
parity eigenstates. 

In the above the argument of Besscl functions is jni—\(l\'>':^)- Using jo,2{~x) = jo,2{x) and 



x) = ^ji{x) we can replace the argument of the Bessel functions to be 



and gather 



together the general structure of the matrix elements for the various polarisation states 



(hyb; ±, m'| exp 



-iq ■ r- 



a ■ x± |conv; I — 0) 



(19) 



and 



(hyb; ±, m'| exp 
1 [2 



-iq ■ r- 



a ■ zjconv; I = 0) 



{5{V = +){f{Jo)^+f{J2)^)}Srn',0 



(20) 



V2f3i V 3 

These are the main equations that set the normalisation for hybrid excitation involving any oper- 
ator. To illustrate how they are applied we return to the specific example of El electromagnetic 
transitions. 

The El operator is ean. (|15|l . thus we simply rescale the above expressions accordingly. Noting 
that |e±|/|x±| = l/\/2 we have effectively 

\mQ WLdJ V TT-^ 

The matrix element for the transition from to spin -singlet hybrid(l++) then follows upon 
multiplying this ratio by the reference form (|19() (for V = +), in the \q\r ^ 1 limit 



(hyb;±,m'|a ■ a?+|conv;/ = 0) = |'5(P 



rn', + 1 



giving finally 



M(conv(0-+) 7 ^ hyb(l++)) ^ ( ^ + ^] \q\J ^S^'.+i I r^dr Rl^,^ir) r i?conv(r) 

\mi 1712/ V 37r-^ J ^ 

which is the form used in ref. |22| . 

This applies immediately to the excitation of the hybrid in jn^ — > afj^ where there is no spin 
flip between the spin singlets tt and ain- Note that this transition requires charged states, the 
neutral modes vanishing in accord with charge conjugation. 



7V/J(77r ^ ain) 

MijTT ^ bi] 

So the ratio of widths becomes 



ei + 62 



ei - 62 



H{r)-n\qH 



b{r)^\qb\ 



2b 
3^ 



2\/3 



IH 



TEi{b 



■ 7r+7) _12 b 
7r+7) TT^ 



19^1 exp 



(_Wh\^\ 



Igbl'^exp 



(21) 
(22) 

(23) 
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where the factor in square brackets includes the g"^ phase-space and a "typical" form-factor taken 
from the case of harmonic-oscillator binding [2^] • These factors model the non-trivial hybrid meson 
mass dependence of the width'* . 

Compare the form of the ratio of El widths, (after removing a factor of 9 due to the different 
charge factors), with the transverse contribution to the elastic charge radius, ean. HlQ(l . In the 
approximations that we have used here, the El transitions to the leading states saturate the 
dipole sum rule. 

Our calculation of the relative strengths of the matrix elements for hybrid and conventional El 
transitions, eans. H21l I22|l . in the flux-tube model, suggests a way of calculating this more directly 
in Lattice QCD. The essential features of the electromagnetic matrix elements are 

(i) initial charged tt^ ; 

(ii) El transition operator with no spin-flip; 

(iii) compute matrix element to = 1+ final states with G-parity = ±1. 

The transition to G-parity -1-1 is to the conventional axial meson, while that to G = — 1 is 
only accessible for the hybrid configuration. If further one calculated (r^) for the tt^ one could 
assess how well the sum rule is saturated by these states, quantify the "penalty" for exciting the 
gluonic modes or hybrids in general and potentially assess the role of such configurations in the tt 
wavefunction. 



2.3 El Rates 

In the Isgur-Paton adiabatic model with a variational harmonic-oscillator solution^ we obtain 
|/f (r)7r/fc(7')7rP ~ 1-0, so the radial moments do not suppress hybrid production. We follow ref. [S] 
and use the standard parameters b = 0.18GeV^,m„ = 0.33GeV so that the prefactor '^-^ ~ 3.8 
and hence there is no hybrid suppression from the flux-tube dynamics. 

Within our variational solution (5^ — 255MeV, /?(, ~ 281MeV, (5-^ = 335MeV, so we see the p = 1 
hybrid state being roughly the same size as the L = 1 conventional state. The main uncertainty 
is the computed size of the tt [23 • Assuming that this hybrid has mass ~ 1.9GeV piinill9j. and 
using the measured width r(6^ 7r+7) — 230 ± 60keV |2] we predict that 

r(a+^ 7r+7) = 2.1±0.9MeV. (24) 

where the error allows for the uncertainty in j23ll24j . We present in Figure.® the dependence 
of the radiative width on the mass of the hybrid. 

The equivalent El process for = 1 is hju — > p7, where the only difference from the S — Q 
case is the addition of L, S Clebsch-Gordan factors coupling the QQ spin and flux-tube angular 
momentum to the total J of the hybrid meson in question. As above, the charge conjugation of 
the initial and final states is the same, thus AG = 0, and the amplitude is proportional to ei -f 62- 
Consider absorption of a positive helicity photon. The hybrid state is constructed as follows 

|J+",TO,/) = ^ (l™'; lm5|Jmj)^{|hyb;-|-, to') - |hyb;-,m')} IS" = 1,7715) 

7715, m' 

and the matrix element becomes 

M{-ip^hjH) = (l+l;lTOp|Jm,7) [ £l±£i j ^(j-) 

\ ran ) 

We find for J = 0, 1, 2 in this El limit, normalising against measured r(/i — > fn(\ 

r(6+^ p^-i) = 2.3 ± O.SMeV. (25) 

^Keeping the full spherical Bessel functions in the radial overlap gives a form-factor very similar to this. Even cal- 
culated true to the model, the form-factors should be considered to be at best rough guides to the mass dependence 
especially in regions where the non-rel approximation is failing 

^see Appendix IdI 
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TO(l++)/MeV 

1400 1600 1800 2000 2200 

Figure 2: El width as a function of 1+"*" hybrid mass. The sohd hne is for = 335MeV. The 
dashed hne is for — 540MeV |23 ■ The shaded grey areas are the uncertainties due to the error 
on the experimental rate used as normahsation. 



where the error reflects the uncertainties in the conventional El strength and and where we 
have taken mn — 1.9GeV. 

We present in TableEl formulae for El radiative matrix elements between conventional and hybrid 
states. Of particular interest is the rate of production of the isovector 1 ^{ttih) at 1.6 GeV. We 
use the El decay of this state to a2 as an explicit example of the use of Table. |21 With a positive 
helicity photon there are three helicity amplitudes, corresponding to mj = —2, —1, 0, 

^(1 - 1; 1 - 1|2 - 2)(10; 1 - 1|1 - 1)M = ^/f-M 

|((10; 1 - 1|2 - 1)(1 + 1; 1 - 1|10) + (1 - 1; 10|2 - 1>(10; 10|10)) = ^^f^t 
^ |((1 - 1; 1 + 1|20)(10; 1 + 1|1 + 1) + (10; 10|20)(1 + 1; 10|1 + 1)) M = ^^M, 
where M ^ \q\ 87 x lO^^GeV^^ Using eqn. 1141 we obtain 

„ , . . Ear, , ^ 1 /, . . ,T 




r£i(7riH 027) = 4^^|9l^(|X_2|' + |M-i|' + |A^o|') 

= -^kf fi + i + -^)(87xlO-3GeV-i)2^ 
m^i„ \2 4 12/ 

so that for a ttih at 1.6 GeV the width is ~ 90keV. Given that 02 exchange is suppressed 
relative to /2 in photoproduction, this is unlikely to be a major production route in jp — > ttihti. 
Pion exchange provides an opportunity via the Ml multipole but this goes beyond our present 
discussion. If calculated via VDM from the ttih — > 7rp rate, where the photon-rho flux-tube is 
excited by a "pion current" we find |2S| r(7ri^f — > -fTr) ~ 50 — 200keV which is similar to that 
found in [EI. 

A state potentially interesting in heavy-flavour decay (see next section) is the axial hybrid Kaon, 
we find that this state has an El width to of 300 — lOOOkeV (assuming mass ~ 2GeV). This 
state could be seen in photoproduction by looking in the /vttttA end-state. 

Note that these El transitions are only possible with charge exchange and so cannot occur between 
fiavourless states. In particular they are absent for cc and bb. Thus, for example, the transitions 
■0(3685) — > 7Xj can receive no contribution from any hybrid component of the V'(3685) wavefunc- 
tion (assuming here that the xjh states are 3> 4 GeV in mass and so do not mix measurably into 
the cc states.) 
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Table 2: Photon-Meson- Hybrid £^1 matrix elements: M = + \^\I^H{r)i should be 
multiplied by the Clebsch-Gordan factor in the second column to give the overall matrix element for 
a positive helicity photon. The numbers quoted in columns three and four are A^/|g| (10~^GeV~^), 
evaluated using the results of Appendix^ except those in brackets which use the /3- values of |24) . 



This illustrates the principal for calculating matrix elements for hybrid excitation by currents. 
We can extend this now to vector and axial currents of arbitrary polarisation and apply to the 
production of hybrid mesons in the weak decays of heavy flavours. This requires a non-relativistic 
decomposition of the currents (Appendix ^ , then we identify the terms linear in a which can 
cause transitions between conventional and first excited hybrids. Finally all one then needs to do 
is to read off the relevant operator, scale by the reference operator as above, and determine the 
relevant matrix elements. 



3 Heavy Flavour Decays to Hybrid Mesons 

In the first half of this paper we have described the method of calculation used to compute ampli- 
tudes for processes involving a hybrid meson, a conventional meson and a current by considering 
the specific example of hybrid meson radiative decay. The results obtained are applicable to the 
forthcoming hybrid photoproduction experiments at Jefferson Lab j26| . 

We will now use the same techniques to calculate the rate of production of hybrids in heavy flavour 
decays. Such a calculation is timely in view of the orders of magnitude increase in statistics on a 
wide range of exclusive decay channels anticipated at present and upgraded B and charm factories. 
We will consider in particular the supposed excess of events in the inclusive decay B J/ipX at 
low J/V'-momentum |13| as this may have an explanation in terms of hybrid Kaon production. 

In Tables (O and we identify the exclusive channels in which hybrids may be clearly observed. 
We begin by outlining the model used to describe the non-leptonic weak decay process. 

3.1 Naive Factorisation Model 

The matrix element for decays B M1M2 will be written in the generic form (see e.g. Ref.( |27p). 

{A'hM2\H,s\B) = ^Vg,bV,>g,0 (26) 

with 

O = ai(^)(Afi|q^(V^-A)6|B)(M2|g2(V^-A)g^|0)+a2(M)(M2|g2(T^-A)fe|B)(Afi|g'i(T/-A)<7i|0). 

(27) 
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The model arises from performing the QCD rcnormaUsation of the weak interaction, supposing that 
one of the mesons is created from the vacuum by a current with the same quantum numbers. Final 
state interactions are ignored, ai, 02 are considered as phenomenological parameters determined 
by fitting to known conventional decay rates and not as the Wilson coefficients they would be in 
the strict theory. 

Decays proceeding by the first term only are labeled "Class 1", an example being — > 7r^Z3'^*)+ 
which in this model has amplitude, 

^VbcVudai{7r-\dA^u\0}{D^*^+\c{V'' - A'')b\B"). (28) 

The pion creation current is conventionally parameterised by {t:^ {q)\dA^u\0) = if-^qp, and we 
compute the hadronic matrix element using our non-relativistic quark-flux-tube model. 

Decays proceeding only by the second term in lf?7|) are labeled "Class IF or "color-suppressed" 
and include the well tested channel _B+ K'^*^^ J /ip with amplitude 

^VbcV,sa2{Jli^\cV^c\Q){K^*^+\-s{V^ - A^)h\B+), (29) 

where the J/ij} current has Lorentz decomposition (J/?/;(g, e)|cVJjc|0) = e*^f^m^. 

Models of this type are more thoroughly discussed in where they are demonstrated to success- 
fully predict weak decay rates of the B-meson to a wide range of conventional exclusive channels. 

We use meson state normalisations as in |28| such that a factor y/AmBmM appears in all hadronic 
matrix elements, and where M is the meson not created from the vacuum. With this normalisation 
the decay widths are written in the following manner, 

r(S°^p-i?+W) = ^-|^|y,,K.nai|' E l-^l' (30) 

^ hclicitics 

where M = {p-\dV^u\0){D'-*^+\c{V'' - A^)&|i?°) and where to be specific we have shown the case 
of the Class I decay to a conventional vector. We will use oi — 1.05, 02 ~ 0.25 for all i?-decays, 
these numbers are compatible with those in and give a good fit to the conventional decays. 

3.2 Application to B decays producing Hybrid mesons 

In this section we extend the established model described above to include hybrids in the end 
state, in much the same was as we extended the conventional radiative decay formalism previously. 
We simply include the additional flux-tube transverse degree of freedom, a, in the currents and 
wavefunctions. As in the radiative case this both modifies conventional decay rates (by changing 
the charge radius or Isgur-Wise function as detailed in Section^]) and allows for hybrid excitation. 

Previous models of hybrid production in this type of process [221 , IHD] have assumed that hybrids 
are created by a color-octet current with some undetermined strength; we differ somewhat in that 
the flux-tube model has had color degrees of freedom removed from the beginning. We have the 
hybrid produced by the same current that produces the conventional states. Rather than color- 
octet and color-singlet factors controlling the production of hybrids and conventional states, we 
have the hybrid production rate determined by excited flux-tube mode in contrast to ground state 
flux-tube mode for conventional states. The relative rate for hybrid production generically is then 
a factor 0{b/m?). 

When considering weak decays of heavy-light systems we use the following form of the quark 
position vector, 

^ ^ /i 2/3ir p2}) 

rQ{d) = R±r ha ■ ' 

mQ(d) niQ + rud 
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This differs slightly from the form used previously, ean. (|12|l . It has the advantage that under 
a finite change in heavy quark mass (i.e. a flavour changing transition) f, which here is the 
interquark separation vector, and dp are unchanged. For details see (31) . 

3.3 B decays to Hybrid and conventional pseudoscalar 

The matrix element for such a decay is 

M{B-^ PH) - ifpq ■ {H\{V - A)\B). 

Examples include tt~D+ (Class I) and B° D° nuH (Class 11)^. Thus to compute the 

branching ratio for n~D^ we need the values of the operators q • A and q ■ V between B and 

Dh wavefunctions. In Appendix 1X1 we compute the terms in the non-relativistic reduction of V^i 
and that are linear in a, and which therefore can hence induce transitions between conventional 
and hybrid mesons. 

We find, for the maximally parity violating current, 




1 



TOfc \q\ 



and for the parity conserving current. 



k1 



2r 

'mo 



2mc 



+ a ■ z A a 



— ( -I + „ igi y 

2mb \ TOc 2mcJ 



q- Ah = e 



'mo 



2mc 



(31) 
(32) 



having approximated q^ — w 0, as appropriate for B decays. 
Spin-singlet Hybrids, L>ij(l±(±)) 

Since {a) ~ for transition to a hybrid where the qq spins are coupled to 5 = 0, the parity 
conserving current ean. (|32|l plays no role here. Non- vanishing contributions come entirely from 
the first term of the parity violating current, ean. 1)31(1 . We thus have a matrix element. 



X(B"^7r-i?+(5-0)) = z/,) 



2b 



{Dh;S = 0|e-*«"^^W— /3i|gla • £ 



-(i+^ + M-] + \q]—(i+M 



2mb 



rrir 2m 



mo 



2771, 



The integrals over a and angles 9, (f) are performed in Appendix IbI and we take the result ean. (|20() . 



(hyb; V = ±, Tt'j exp 



-iq ■ r- 



a ■ zjconv; I — 0) 



V2/3i V 3 



^recent theoretical work |32| suggests that naive factorisation can be a poor approximation for processes 
1?'^+ light hadrons 
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where in the current example the argument of the spherical Bessel functions is \q\r . Thus 

we have the matrix elements, 

A^(B° ^ 7r-D+(l-(-))) = 0, (33) 

[ 2m6 \ mc 2mc J nid \ 2nic J ) 

That the l^(^) amplitude is zero was to be expected as one cannot maximally violate parity and 
conserve angular momentum in a process 0^ 1^0^ in any partial wave, whereas the 
amplitude is non-zero as 0" 1^0^ in a P-wave respects the symmetries. 

Spin-triplet Hybrids, Dh{J^^"^^) 

To excite spin-triplet hybrids we require operators linear in S and a, which feature in both q ■ A 
and q ■ V, eans. (|31l32|l . 

For simplicity in presentation we shall define p± = ^1 ± ^ -f 2^^) ■ We can decompose 
a • z Aa = ^(cr+a • x_ — (T_a • x+) and a ■ a = ^{a+a ■ x_ + (7_a • x+) + (JzO. ■ where cr±, cr^ are 
normalised such that {S = l,ms\VAS = 0) = ["^"f J"of • 

The = 1 of the quarks and L = 1 of the flux tube combine together to give the J of the 
hybrid meson, with the appropriate Clebsch-Gordan coefficient (Imi; lms\Jmj) . Performing the 
integrations over a and angles 0, as in Appendix Ogives for the parity violating current. 



Mq.v ^ -fnV^m^nwJ ^^\q\p- ^ (lmL;lms\Jmj) X 



37r3 

(35) 



S{V = +) {(jo) - l{j2)) {SmL.-lSms. + l " SmL. + l^ms-l) 
+ S(V = -)f^(jl) {SjnL,-lS,ns, + l + SmL, + lSras,-l) 

The pattern of amplitudes for various J^'~' follows from the combination of Clebsch-Gordan coef- 
ficients, 

SiV = +)(Oo) - 1{J2)) [(1 - 1; 1 + 1| JO) - (1 + 1; 1 - 1| JO)] 

+S{V = -)|(ji) [(1 - 1; 1 + 1| JO) + (1 + 1; 1 - 1| JO)] 
Thus we find, for the maximally parity violating current, 

o+(-) 

2+(-) 

An exactly analogous process for the parity conserving current gives, 

Mg.A 
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( j: ) . (a { 3> ^ ^ (- i:) - { i )) 

(37) 

As in the spin-singlet case we can understand the zero values as originating in the need for angular 
momentum conservation while applying the appropriate parity selection rule. 

Numerical Estimates 

The model parameters and variational wavefunctions used to evaluate the radial matrix elements 
(DhUlIB) are described in Appendix IdI We obtain the following results for the branching frac- 
tions of the meson to the exclusive channel D~^{J^'^'-^^)n~ , presented in Table.©- 



uidh (GeV) 

Oo) 
(.7.) 


2.7 

0.62 
0.24 
0.06 


3.0 

0.65 
0.22 
0.05 


3.3 

0.67 
0.21 
0.04 




l + ( + ) 


6.1 


4.5 


3.0 


xlQ- 


-4 


l-(-) 















0-(+) 


3.1 


2.5 


1.9 


xlO" 


-7 




1.3 


1.3 


1.2 


xlO" 


-6 


2-(+) 


1.6 


1.3 


0.9 


xlO" 


-7 


o+(-) 


3.0 


2.3 


1.6 


xlO" 


-4 


l-(+) 


4.9 


3.6 


2.5 


xlO" 


-6 


2+(-) 


3.3 


2.3 


1.5 


xlO" 


-4 



Table 3: &.r.(BO ^ n- D+iJ^'^^^)) 

Whilst these are the decays in which the naive factorisation approximation is most likely to be 
correct they are not ideal for hybrid hunting: the D-mesons are not eigenstates of C and as such 
cannot have "exotic" quantum numbers, the smoking-gun signature for a hybrid. However, the 
fact that for a given J^^'-^^ the hybrid and conventional states have qq coupled to S' = 0(1) and 
S — 1(0) respectively can lead to selection guides. This has been noted already for hadronic 
processes jHIESl- We shall see that there are further examples in B-decays. 

The Class II decay 5° — > D'^{nn)H, by contrast, can produce exotic quantum numbered hybrids 
and is not suppressed by any small CKM matrix elements. Unfortunately there are problems both 
with factorisation j32| and with our model formulation when applied to this decay. Our model 
has form- factors obtained from the simple wavefunction overlaps and although these 

are reasonable at small |g| where everything is non-relativistic they are at best a qualitative guide 
as the phase-space rises. Such subtleties are discussed in where a better fit to the spectrum 
of semi-leptonic B decays is found using a power law form-factor in contrast to the "polynomial 
times exponential" form that we have used. To reduce the uncertainty we calculate the ratio of 
hybrid production and conventional (known) rates, for which much of the form factor dependence 
cancels out. Thus for example, if we propose that there is a 1 hybrid at 1600 MeV we predict 
that it will have a branching fraction a little smaller than that of i?° — > D^{p,uj) i.e. O(10~^). 
(0, 2)^ exotics have potentially even larger rates provided that they have masses somewhat below 
3 GeV. 

In addition to B-meson decays, this formalism can equally well be applied to D and Dg decays, 
though the factorised model may not be as accurate here. There is some possibility of producing 
exotic hybrids in the Cabibbo suppressed decay Tr~^ {nn)Jj if the hybrid is light enough, 

similarly in the channel Dg 7r+ (ss)ij. The numerical branching ratios for these processes are 
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strongly dependent on the hybrid mass as they are so close to kinematic threshold. The candidate 
1 ^ state at 1600 MeV has a branching ratio in this decay of 0(10^^) principally because it is 
produced in a P-wave with a very small \q\. 

3.4 B decays to Hybrid and conventional vector 

The matrix element — > VH) = my/ye* — A)^\B), where is the vector meson 

polarisation. We thus have a set of amplitudes depending upon the helicity of the vector. Examples 
of this type of decay include p^D^*^^ (Class I) and the well tested decay J/^pK''*^^ 

(Class II). 

As discussed in Section^ the |A^p for the weak transition B ipKnil^) is expected to have 
strength ~ 13% relative to its "conventional" counterpart B ipK{l^). Empirically 5+ ^ 
■0iv:(l+)(128O) is the single largest mode in B+ ipX with b.r. = (1.8 ± 0.5) x lO'^ while 
5+ -0i4r(l+)(14OO) < 0.5 X 10""^. These rates involve both parity conserving (vector) and 
violating (axial) contributions and their relative strengths depend on the mixing between the 
•^Pi and ^Pi basis states. These rates would lead one to expect an order of magnitude b.r. for 
5+ — > ^PKh{^^) > 10"'*. This has been reported by us in ref. Here we illustrate the 

calculation. 

Spin-singlet Hybrids 

As previously we select the terms in the hybrid currents (Appendix independent of a, for which 
we find terms in both and A'^. We compute separately the amplitudes for longitudinal and 
transverse J/ip'. 

Longitudinal J/tp and Kh{1^) 

For a longitudinal J/ip, = (I*?!, 0_l, P^)/™,/, and so 

Ml = U{KhMV^-E^V^\B) 
since only Ah± are non-zero for spin-singiets. Using the overlaps in ean. (|2U|l gives 

{ md ( ^ 2ms) ^'^^^ ^ 2mf, (2ms ^ ''^ (m^ ^ )) ^^^'^^ 02))| (38) 
Thus only the positive parity state is produced; there is no longitudinally polarised if// (1^^^'). 

Transverse J/t/i and _K'//(1*^*)) 

For transverse polarisation e*^ = (0, e*) and e* — t(1j T*, 0)/\/2. Both vector and axial can now 
contribute. The vector current leads to a matrix element 

Ml = z/^mv,VW^yj^ (^^ + 1) <5,„, ±1 |<5(P = +) (^(jo) - ^(^2)) ± 5{V = -)|(ji)} , 

(39) 

and the axial current to 

Mi = ±z/^m^V4^^I^yj^^<5„-,±i l^SiV = +) (^(jo) - \{J2)^ ± |'5(P = -)(ji)| 

(40) 
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Hence the and transitions combine to give the total matrix element for transition to !+(+) 

(41) 

The transition to 1 ^ ^ has the same structure apart from the partial wave contributing as ji, 
and reads 



(42) 

We see that for each of these transitions, the vector current dominates. Hence we may anticipate 
that the transition to Kh{1~^) will be relatively large because the dominant {Vfj) contributes in 
S-wave; by contrast Kh{^ ) receives its S'-wave from the |g|/m;,-suppressed (A^) while the vector 
current contributes to P-waves. Explicit calculation, below, confirms this. 



Longitudinal to transverse ratio . The 1~ hybrid is produced in transverse polarisation 
only. For the 1+ hybrid both transverse and longitudinal polarisations are possible. Compare 
the longitudinal matrix element, eon. H38|l with the transverse ean. H41|l and remove common factors 
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2m,, 



(jo) 



\{J2) 



1 



\ms 
2m, 



(A) 



Before taking the ratio, consider the actual values of the parameters. For a 1.9 GeV Kaon hybrid 
we have \q\ = 0.83GeV, thus \q\/2ms ~ 0.8 is negligible next to (1 + mt/ms) ~ 10.5. The J/tp 
is moving non-relativistically such that for this \q\, w is good. Using our variational 
wavefunctions f Appendix IDl) we find (jo) ~ 0.58, (ji) « 0.23 and (j2) ~ 0.06, so we can safely 
neglect (j2). Thus we find. 



Ml 



Mn 



1 + ]^J^b |g| 



1 



TT m^ md V 1 + mii/m 



m^/2ms\ (ji) 
Oo) 



3.4, 



which leads to a longitudinal width fraction. 



%Tl = 100 



100 1 



Ml 



Ml 



90%. 



This number should not be taken too seriously however. It has proven difficult to accommodate 
the measured value of the longitudinal width fraction of 5+ J/4>K*^ within factorised models 
and this failure may signal the limitations of this simplistic model. 
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Spin-triplet Hybrids 

Considering the terms in eans. l49l51l53l5H)l linear in a yields matrix elements, 
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Ml I 
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^^jUn) } 1(43) 
71 



Ml I M = /v.V4to^to^\/^ X 
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^ 2mb \TOs 2msE^ 



(Ji 



(44) 



l+(-) 
2' 



M±[ ] = iUv'4:mBmK\l —— x 



26 1 
±1 

2m., y I 1 ) 2mf, V m 



2m< 



2(jo) + i02) 
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M- 



l-(+) 

_(+) ) = UV^mBm^' 



'2b 3 



2mb ^rUs 2ms) l 1 



(46) 



We can use these matrix elements in the width formula 

^2 

IGtt m 



r(B+^J/V.if+(J^(^))) = ^J-|HeK.na.p 5: = 
to compute branching fractions. 



Numerical Estimates 



We explicitly evaluate branching fractions and longitudinal rate fractions for the exclusive channel 
Jf-ipK^ using the parameters and wavefunctions given in Appendix (DJ the results being 
presented in Table. 

While fine details of the model may be questioned, the O(10~^) branching ratio to the hybrids 
with positive parity appears robust and accessible to experiment. It is intriguing therefore that 
there is an unexplained enhancement at low q^, corresponding to high mass K systems, of this 
magnitude [O] . 

While suggestive, it would be premature to claim this as evidence for hybrid production. Radial 
excitations of the K(l^) are expected in this region, and in the ISGW model, extended to 
exclusive hadronic decays and assuming standard factorisation arguments |27j . we find these to 
have b.r. ^ 10^"*, though slightly less than the hybrid. 

Other conventional strange mesons in this mass range are likely to be suppressed due to their high 
angular momenta, which give powerful orthogonality suppressions at small \q\. It is the S'-wave 
character of the hybrid and axial production that drives their significant production rates. 
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Table 4: 6.r.(S+ ^ J/iPK+{jP<-'^^)). The subscripts are the longitudinal rate fractions. 

The channel i?" — > D'^*{nn)'jj, whilst suffering the same theoretical problems as i?" — > D'~'nn^ 
considered earlier, is not in principle suppressed by any small numbers. CLEO II has observed 
events B'^ Z?"*7r+7r^7r+7r~ ^3] which is a possible end state for light-quark hybrids decaying 
via af?:^ for example. According to the flux-tube breaking model of hybrid hadronic decays |S] 
isovector 2^ , 1^ , 0+~, 1 and isoscalar 1 ^, 1++ hybrids have large branching ratios to ain. 

4 Discussion 

The results of this paper suggest that the flux tube can be excited without undue penalty. The 
physics assumes that a flux tube is indeed formed and if subsequent lattice studies should confirm 
that the 1 fm length scale of confinement is intermediate between that of perturbative or bag-like 
gluonic fields and a fully formed flux-tube, then some of our results may need to be reassessed and 
more mature modelling developed. However, this approach seems likely to exhibit features that 
will survive, at least in part, in a more mature description. 

There is an obvious question as to the reliability of our non-relativistic treatment for light flavours. 
We suggest however that the physics of the excitation is probably more general than the specific 
modelling herein. For very massive quarks, where the non-relativistic treatment may be justified, 
the cm. of the system tends to lie on interquark axis r = rg — rg. For \q\ < mg the recoil of 
the quark is small and the cm. tends to remain near this axis. As a consequence, we find that 
excitation of the tube is suppressed. As the quark masses are reduced, the cm. of the system 
can have increasingly large excursions from the r axis, and the tube is more easily excited. There 
comes a point where the quarks are light and the non-relativistic approximations are suspect; 
however, the physical picture that the flux-tube excitation is easier for this situation than for the 
heavy quark case is a physically reasonable extrapolation. So although the actual numbers may 
be debatable, the probability for excitation of the flux-tube seems likely to be at least as big as 
was calculated for massive quarks where the non-relativistic approximations still apply. The J^'~' 
patterns of the results, and the conclusions about charge exchange processes being required for 
the excitation of hybrids in the El multipole, also seem to be robust. 

Our results suggest a way of assessing the flux-tube excitation in a lattice QCD approach. For 
an El transition with no spin-flip, a tt can be excited to axial mesons. The electrically charged 
hybrid {a^j^) and conventional (bf) axial states thus excited have opposite G-parities. Thus the 
relative strength of the matrix elements for electromagnetic El transition to G-parity ±1 gives a 
measure of the penalty for exciting the flux-tube. 

It will be interesting to see if lattice QCD simultaneously can describe the magnitude of (r^)^, 
the El amplitudes for exciting G = ±1 axial mesons, and test the extent to which the electric 
dipole sum rule relating these is saturated. In our approximations we see that it is saturated by 
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the ^Pi conventional states and their hybrid axial counterparts. If these results are verified, then 
it may be possible to extrapolate to other multipoles, relating the modifications of various static 
properties by the flux-tube degrees of freedom to the excitation of hybrid states of various J^*^ . 

The hybrid configurations with non-exotic J^'-^ will tend to mix into the wavefunctions of the 
conventional mesons. Unless specific observables (such as ^) have anomalous values, it will be 
hard to make a convincing case for such states. If the enhancement in _B ^ can be shown 
to be due to a Ka state, this would be very interesting but would not of itself be proof of a 
hybrid. Further measurements, of polarisation or ^ would be needed to distinguish it from a 
radial excitation of a conventional axial. The main interest is seeking rates for states with exotic 
J^'-'' . Table. 121) and eon. 1)25(1 predict healthy couplings to jp, which implies that photoexcitation 
of 0^ , 2^ off a p is feasible. We note also that diffractive production of 2^ with photon beams 
is to be expected, hence 2"' exotic states should be searched for in charge exchange and also in 
diffractive scattering e.g. at either Jefferson Lab or HERA. This is discussed further in 125,1 . 

In the decays of heavy flavours B — > DhX hides the exotic quantum numbers in the flavour of the 
Dh- However, the "spin inversion" between conventional and hybrid vector mesons (the hybrid 
being a quark spin singlet in contrast to the conventional triplet) leads to an interesting selection 
rule: 7r^L)j!^(I^(^') = 0. Hence observation of a vector D* in other processes, that is absent 

in B decay, could be a signature for Dh{1^). 

The production of an axial strange hybrid in i? — + ipKn is predicted to have b.r. ~ 0(10^^) so long 
as its mass < 2.IGeV. There is a tantalising unexplained enhancement in the data that may be 
compatible with this and merits further investigation. This process superficially has an analogue in 
B D*^nnH, which opens up the available phase space for nfi to ~ 3GeV enabling light flavoured 
hybrid states with exotic quantum numbers to be produced. Unlike the previous case, however, 
there is the possibility that rescattering effects and the failure of factorisation could contaminate 
our analysis of this process. Modulo this caveat, we have b.r.{B D*^(nn)H) ~ 0(10"^) for 
7ri(I600), if this is the exotic hybrid with J^*-^ = f ^; (0, 2)^ exotics have potentially even larger 
rates provided they have masses below 3 GeV. 

The exotic f ^ resonance 7ri(1600) could in principle be looked for in D decays. We find a 
branching ratio for D 7r7ri(I600) that is ^ 10^^ times that for D — > np. Folding in the 
predicted PSIISS] &.r(7ri(I600) np) ~ 20% gives a combined branching ratio to mrp via this 
exotic of ~ 10~*. This would be a severe challenge even for high statistics studies that may become 
available at CLEO-c or GSI. 

This exotic should have 1=0 partners, if it is indeed a hybrid and not some di-meson effect. The 
{&)h state can be produced in leading order Ds ttssh if its mass is below 1.9 GeV. If the nonet 
is not ideal then there is the hope of some fi content in each of the 1=0 states, one hopefully 
light enough to be accessible. Similar to the above, we find a branching ratio that is ~ I0~* times 
that for Ds TTcf). Here again, one is at best at the limits of detection. Another potential source 
of these exotic hybrids is in the decay B^ J /^pssH which is the s-quark spectator analogue of 
B J/ipKn- We find a branching fraction for a I ^ ss hybrid at 1.9 GeV to be ~ 6 x 10^*^ 
which is far lower than current statistics can observe. 

The general formalism developed here can be applied to any current induced transition. Examples 
include gluon emissions (as in cc cascades), diffractive excitation or tt emission where the pion 
is treated as an effective 75 current. The latter has a long and successful history in describing 
conventional hadron decays |36ll37j and can now be applied analogously to the production or 
decays of hybrids involving tt. This is described in ref. [221 
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Appendices 



A Hybrid transition operator for V/j, and 

We perform the non-relativistic reduction of the vector and axial currents allowing for flavour 
changing. We present the particular case of B{bd) D{cd) by quark level b ^ c. In the following 
mo = md + rric and ms = + m^, which are indeed the meson masses in the extreme non- 
relativistic limit. 



For the zeroth (time) component of the vector current we have 



4m6mc 4 \mc rrib 



(47) 



Considering only the terms linear in a obtained from utilising the effect of the momentum operator 
on flux-tube ground state wavefunctions, 



P\Xo) 



= -ti / — /3i a|xo) 

a— cmpt V TT 



(48) 



and expanding the plane wave to leading order in g • a we have the effective transition operator to 
the first hybrid excitation 



' ■ n. \ 



/3i { iq- a 



+ 



-a ■ qA a 



(49) 



V 



For the spatial vector current we have 



V{pc,Pb;x) = c{x)jb{x] 



J{Pc-Pb)- 



Pb . p, 



+ 



2mb 2m, 

and the relevant transition to first excited hybrid becomes 



-ia A 



Pb Pc 



2mb 2mc 



+ .. 



VH = e 



-iq-r- 



2mb Vm, 



rricmD 



2mb \m, 



(50) 



'2b ^ \ .TT /nib ,\-, . r ^ tt /nib -> r 

rPi 1 ^^t; — I 1- 1 ) a — ^ [q ■ a) q+ - — ( 1 \ a Aa-\ [q ■ a) a A 



nicniD 
(51) 



A°{p,,pb;x) = c{x)j''^%{x) e^ff^-f")-^ ^ a ■ 



Pb _^ Pc 

2mb 2mc 



and the relevant hybrid transition becomes 



TT / nib 



2mb Vm, 



1 ] a • a — i- 



rricniD 



-(9 



• a) a • a| 



(52) 



(53) 
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A 



A[p^,p,,x) = -c{xWb[x) ^ e^(---)- + ^'^^'^ ' + ^^^'^ ' zfMl 

(54) 



and the relevant hybrid transition becomes 



Ah '"D ^ —fj^ <i{q-a) a [ + i- [(a ■ a)q + {a ■ q)a\ + q A a 

(55)' 



B How to calculate transition to hybrid 

We consider the generic 

M = (hyb; ±, m'|acf |conv; l,m) = J (fr j (fa TCO^^fC 

where O-^ci = a ■ Xi and the ± is the hybrid state flux-tube polarisation. In this specific example 
we shall choose = x ~ iy and calculate the matrix element 

(Xi,±|a • x_|xo) = -(Xi>±l(ai + «a2)2?i+* - (ai - ja2)2?-l* |xo) 
where we have introduced the T) rotation functions shown explicitly for j = 1 in Table ^ . 

+ 1 / i(l + cos6l) -^e-*"*sin6l ie-2«'^(l - cos^ 
^e*'^sin6' cos 6* -^e"'"^ sin 61 

-1 yie2^'^(l-cos6') ;^e''^sin6' i(l + cos6') 

Tables: (0,^,-0) 
The explicit expression for the matrix element is 

d?api{aiT^a2) ((ai + ia2)I?i'|* - (ai - laa)!?^'!*) |xo(ai)nXo(a2)|' 
where the xo a-re as in jSJl- 

Rewrite ai ± ia2 = a* and the integral becomes 

-P, j d2aaT(a+I?it-«"2'-'-)IXo(ai)nXo(a2)p. 
Now use a^a~ = a\+ a\ and note that a* a* vanishes under integration. This brings us to 
-/3i / d2a(a2+a2)(^+p(i|*_^-pW*)|^^(„^)|2|^^(^^)|2 



Then using the integral 

?2 
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we finally obtain the essential angular decompositions as follows: 



(Xi,±|5-f-|xo) = (56) 
(Xi,±|a-f+|xo) = (57) 
(Xi,±|a-i|xo) = --^(<5+I?«*-ri?(L)*) (58) 



To proceed from ea ll8l first expand the exponential in terms of partial wave angular states 



'YliL *^(2i + 1)2^00 contract together the three V functions 



— (LO; 1 -1|1to') ((5+(L0; 1 +1|1 +1) - (5-(L0; 1 -1|1 -1)) . 
and integrate / d£l^ which gives for the matrix element 

and finally 

'i^T!'^"'"' [(/Oo).-i/02).) ('5+-<5-)-^^/0•l),(<5+ + r) 

where we now only need to calculate the radial expectation values of the spherical Bessel functions. 



C "Dipole" form in the adiabatic flux-tube model 

In the equal mass case we have for the quark positions, 



,1- br 2 1^ 
2 nm V A* + 1 '^-^ p 

p.odd 

Specialising to p = 1^ we find the momentum of the quarks at lowest non-trivial order in — , 



PQ^d - ±Pr - V2{N+l)ps (60) 
where pp,ps are the momenta conjugate to r, a. 

In the usual non-relativistic quark model the matrix element of the lowest order electric dipole 
operator ~ ^ can be transformed using the commutator p = im[H, r\ into an explicit dipole form 
^ m|g|e*- r. This commutation relation is valid provided the system Hamiltonian can be written 

2 

in the form H = + V{r), any other dependence of p causing deviations from this behaviour. 
We find that this transformation is justified in the adiabatic approximation to the flux-tube model 
also: 

Minimal coupling of the photon flcld to a quark at rq moving with pq leads to a convection 
current operator ^ pq ■ A{rq). For a plane- wave photon field at lowest order in q.rq we get an 
El operator ~ We consider the matrix element of this between a ground state meson and 

^we are implicitly making an adiabatic approximation so that the phonon modes do not mix, for corrections to 
such an approximation see 1391 
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a. p = 1 hybrid state with one excited phonon. Only the p = 1 term in the sum in pq can excite 
the flux tube to give a non-zero overlap with the hybrid state, the ps^ acting on the ground-state 
flux-tube wavefunction xoi (^i)xoi (^i) to give a factor i-^^ai. The full matrix element calculated 
this way comes out to be proportional to 



r^drRnRo- 



(61) 



If instead we use the relation pq — imqlH, rg] to transform the operator and expand the commu- 
tator to get a factor {£^(fl-) — £-(0)) ~ (exact in the non-rel limit), we can use the p — 1 term in 
the sum in fg to excite the flux tube. This leads to a matrix element proportional to 



k1 



r^drRnrRo 



(62) 



which appears to be different to what we had previously. If however we assume that the hybrid 
energy is higher than the conventional meson energy by the additional energy in the flux tube, 
(^) only, which seems reasonable in this adiabatic picture where the quarks do not move, then 
\q\^ = 7r(i)-^ ^ 1 * and the two matrix elements approximately agree. 

For example in the adiabatic potential model described in Appendix ||D)| with / — > oo we have for 
the charmonium sector. 



^191 



3143MeV 
4202MeV 
1059MeV 

1157MeV 



(63) 
(64) 
(65) 

(66) 



So we see good agreement. 



D Hybrid Hamiltonian 

We follow the formulation of the flux-tube model described in [3S], [31), 001 ■ In the appendix 
of |38j the authors show explicitly that the spatial wavefunction of a hybrid state with phonon 
occupation {np_|_ , np_ } can be written 



RnLNAir)^^^V'^lli4>, e, -0) 

where the quantum numbers are; n, radial; {L,mL), angular momentum; N — X)p('^p+ + '^p-); 
A = ^p(np+ — np_). We are only interested in the lightest hybrids which have one phonon excited 
in the p = 1 mode and hence N = 1 and A = ±1. In the adiabatic approximation we get a radial 
Hamiltonian (acting on rR{r)) for these hybrids, 

whereas for conventional mesons (no phonons) we would get, 

i/" - — — + ^(-^ + ^) + 5r - - + c 
2/i dr"^ 2/ir2 r 



^(H\l/r\lS){lS\r\H) = Y:„{H\l/r\n){n\r\H)Sn,is ~ 1 if {H\l/r\n ^ 15) Ri or {H\r\n ^ 15) Ri 0. In a true 
HO system (i.e with just one /3 for all states and <5 = 1 for the hybrid) all the matrix elements of r vanish except 
the lowest and the approximation is exact. The deviation from this case is presumed to be not too large. 
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The modified angular momentum barrier in the hybrid case has its origin in the A = ±1 carried 
by the phonon in the tube, br is the mass energy of the string, f is the excitation energy of 
the string in the p — 1 mode, the additional factor multiplying this is put in by hand and is 
designed to model the fact that at short distances we do not expect "stringy" configurations to 
dominate in QCD. The remaining potential terms — + c are introduced by hand, the first of 
which represents one-gluon-exchange dominance at short distances and the second is required to 
describe the observed meson spectrum. 

The parameters, mq^b, k, c are chosen to reproduce approximately the observed conventional 
meson spectrum (up to spin-dependant splittings). We use the following set of values; 

b = 0.18GeV^ / = 1, c = -O.TGeV, m„,d = 0.33GeV, = 0.55GeV,mc = 1.77GeV, = 
5.17GeV. K is allowed to run in a reasonable way so that for light mesons (nn, ns, ss) n — 1.07, 
for heavy-light mesons (nc, sc, nb, sb) k = 0.67 and for heavy-heavy mesons (cc, cb, bb) k—0.52. 

We solve the Schrodinger equation variationally using a Harmonic Oscillator basis 

- ^ r(n+?'1l/2) ^''^''''-''^-^'''(^'^'^^"P-^'^'/'- 

For conventional states L' here is just the angular momentum quantum number L. For the hybrid 
Hamiltonian the modified angular momentum barrier is cancelled if L' is chosen so it satisfies 
L'(L' -t- 1) = L(L + 1)- A2, for L = 1, A = ±1 this means L' = S ^ 0.62. 





L = 







L = 


1 




Hybrid 








2 / 


^^5/2 










7ri/4 Y 


\J r(3/2- 






fiis 


Mis 


Pip 




Mip 


Ph 


Mh 


nn 


0.334 


0.672 


0.280 




1.296 


0.260 


1.825 


ns 


0.370 


0.780 


0.306 




1.386 


0.275 


1.961 


ss 


0.426 


0.857 


0.342 




1.448 


0.311 


2.024 


nc 


0.389 


2.036 


0.331 




2.544 


0.291 


3.128 


sc 


0.469 


2.092 


0.387 




2.571 


0.335 


3.195 


cc 


0.639 


3.129 


0.509 




3.539 


0.426 


4.235 


nb 


0.408 


5.392 








0.302 


6.449 


sb 


0.510 


5.455 








0.355 


6.563 


cb 


0.798 


6.418 








0.486 


7.575 


bb 


1.239 


9.522 













Table 6: /3 values/ state masses in GeV 



The Hamiltonians , are found to be diagonal in this basis to a very good approximation if 
the P values listed in Table© are used. 

Merlin lU' and Merlin & Paton [3^ consider non-adiabatic corrections to this Hamiltonian, the 
values quoted in Table JHl are actually obtained using this modified Hamiltonian, although the 
differences in P with respect to using are usually small. 
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